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Abstract 

The magnetic field of the permanent mag- 
net electron cyclotron resonance (ECR) ion 
source SWISSCASE located at the Univer- 
sity of Bern has been numerically simu- 
lated and experimentally investigated. For 
the first time the magnetized volume quali- 
fied for electron cyclotron resonance at 10.88 
GHz and 388.6 mT has been analyzed in 
highly detailed 3D simulations with unprece- 
dented resolution. The observed pattern 
of carbon coatings on the source correlates 
strongly with the electron and ion distribu- 
tion in the ECR plasma of SWISSCASE. Un- 
der certain plasma conditions the ion distri- 
bution is tightly bound to the electron dis- 
tribution and can considerably simplify the 
numerical calculations in ECR related appli- 
cations such as ECR ion engines and ECR 
ion implanters. 

1 Indroduction 

The reliable and stable operation of an Electron 
Cyclotron Resonance (ECR) ion source in produc- 
ing high intensity beams of highly charged ions 
depends on the successful design of the magnetic 
confinement. This confinement enables the ECR 
process and is responsible for the confinement of 
high energy electrons for the production of highly 
charged ions. The resonance criteria for the ECR 
process is given by Eq. ([T]) [21 [4] . 



Bres is the magnetic field fulfilling the resonance 
condition, uj is the angular frequency of the applied 
electric field, rrie is the electron mass and e is the 
charge of the electron. There is very little that can 
be done to change the magnetic field geometry after 
production of the permanent magnets. It is there- 
fore very important to perform reliable and precise 
numerical simulations of the magnetic field distri- 
bution before the production process to guarantee 
a successful ECR function. 

In this paper we present the simulation of 
the newly developed and operational SWISS- 
CASE (Solar Wind Ion Source Simulator for 
the CAlibration of Space Experiments) ECR ion 
source (Figure [l] and |2| and demonstrate the sim- 
ulation precision that is necessary for the reliable 
localization of the ECR process. In SWISSCASE 
the ring shaped magnets establish the axial confine- 
ment featuring rotational symmetry with respect to 
the beam axis. The Halbach hexapole [13^ (Figure 
[2| establishes the radial electron confinement. 

The combination of the ring shaped magnets and 
the hexapole requires a full 3D simulation, rather 
than two 2D simulations, which would require far 
less computation power. Therefore, a three dimen- 
sional simulation has been performed covering the 
full spatial extent of the magnetic confinement. 

In addition to numerically determining the mag- 
netic field distribution, the simulation also allows 
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Figure 1: A cross section of the relevant compo- 
nents of the operational ECR ion source SWISS- 
CASE at the University of Bern. Arrows indicate 
the direction of magnetization in the ring shaped 
permanent magnets. Dimensions in mm. 



calculating electron trajectories. In this paper we 
show that the resulting electron current density 
and electron charge density coincide with experi- 
mentally observed surface coating and sputtering 
shapes. Further calculations show the ions are 
bound to the same location as the high energy 
electron population. This explains the coincidence 
of the simulated electron distribution and the ob- 
served ion surface coating patterns. 

In Section [2] we present the calculation of criti- 
cal plasma parameters to justify the results of the 
numerical simulation. In Section |3] the results of 
the magnetic field simulation are presented. Sec- 
tion m details the simulated electron distribution 
followed by the experimentally observed ion coat- 
ing patterns described in Section [5] In Section [6] 
the simulation results are compared to our obser- 
vations. 
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Figure 2: A section cut through the Halbach 
hexapole magnet of SWISSCASE reveals the char- 
acteristic alternating orientation of 24 permanent 
magnet sectors. Bold arrows indicate the resulting 
main poles inside the hexapole. Dimensions in mm. 

energy electrons as well as low energy ions inside 
the ECR plasma. It is assumed the plasma is dom- 
inated by a cold electron population [8]. It is fur- 
ther assumed the plasma frequency to be the cut 
off frequency for the incident microwaves inhibiting 
further microwave plasma heating and limiting the 
cold plasma electron density [TJ [2] . This leads to 
a stable equilibrium where the cold plasma density 
is bound to the microwave frequency. Hence the 
cold plasma density can be estimated by equating 
the microwave frequency and the plasma frequency: 
Eq. ([2| [2]. Ah numerical values are for SWISS- 
CASE. 



fpe = fMW = 10.88 GHz 



(2) 



fMW is the incident microwave frequency and 
fpe is the plasma frequency. From the plasma fre- 
quency one gets the density as [U [2] : 



Knowledge of the plasma parameters is required for 
the calculation of mean free path lengths for high 



, = 1.46* 10^^ 



(3) 
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LOpe is the angular plasma frequency, eo the vac- 
uum permittivity, rrie the electron mass and e its 
charge. 

Because of the acceleration of electrons in the 
ECR zone, cold electrons are heated to substan- 
tial energies, forming the hot electron population, 
which is responsible for the ionization of atoms and 
ions. 

Lacking our own Bremsstrahlung measurements, 
which allows the determination of the plasma elec- 
tron temperature, we assume a temperature of the 
hot electron population of 25 keV^ which has 
been measured by R. Friedlin et al. [5] who ob- 
tained their measurements by operating an ECR 
discharge at 7.25 GHz. With a value of 25 keV 
we obtain the mean free path lengths of the hot 
electron population according to Table [l] Assum- 
ing the same hot electron temperature for SWISS- 
CASE, operated at 10.88 GHz, as for the 7.25 GHz 
source provides a conservative estimate for the hot 
electron energy and hence for the collision cross sec- 
tion. 



Parameter 


Value 


Term 


^Spitzer 


2.05 * 10-2"^ 




^Spitzer 


3.33* 10^ m 





Table 1: (ispitzer- collision cross section for Spitzer 
collisions, Ink: Coulomb logarithm, Ue'- plasma 
electron density, Ej^i^: kinetic energy of the elec- 
trons. 

Different formulas for the calculation of mean 
free path lengths according to different collision 
mechanisms can be found in the literature (Table 



Coll. mech. 


A 


Ref. 


Spitzer 


3.33* 10^ 


m 


Geller [2^ 


Single 90° 


5.90 * 10^ 


m 


Chen [1^ 


Semi empirical 


3.54 * 10^ 


m 


Huba [3 



Table 2: Mean free path lengths for different colli- 
sion mechanisms. 

For the source operation with CO2 gas as will 
be discussed in next sections, the following calcu- 
lation shows the collision mechanism between elec- 
trons and neutrals lead to far shorter mean free 
path lengths than predicted by electron ion col- 
lisions. The neutral population consists of CO2, 



CO, O2 and O3 molecules and of O and C atoms. 
We simplify the neutral population to a combined 
total density of 2 * nco2- It is further assumed 
this hypothetical population to feature an effective 
cross section given by a disc with a radius equal 
to half the length of a CO2 molecule (116.3 pm). 
This results in an effective cross section of — 
4.25 * 10~^^ w? . The mean free path can be calcu- 



lated by Eq. Q: 



A 



CO2 



2 nco2 cr'(jo2 



(4) 



Aco2 is the mean free path in the CO2 neu- 
tral gas background and nco2 its density. Based 
on fluid and thermodynamical calculations and as- 
suming the neutral gas temperature to be 300 
the neutral gas density in the plasma chamber is 
4.805 * 10^^ 1/m^ resulting in AcOs = 4.90 m. 
Hence the limiting mechanism for the mean free 
path is given by electron collisions with neutrals 
rather than electron collisions with electrons or 
ions. 

Another effect that must be considered is radial 
electron loss across the magnetic field lines by col- 
lision induced diffusion. A simplified expression for 
the diffusion process of the hot electrons is given 
by Eq.[5]ffl. 



Fg = — /i He E — De VUe 



(5) 



Fg is the loss rate of the electrons per time and 
area unit, ji is the electron mobility, Ue the electron 
density, E the electric field and the electron 
diffusion coefficient. is given by Eq. [g] [2 . 



(6) 



is the Larmor radius of the hot electrons and 
Vn-e the dominating collision frequency between 
electrons and neutrals as shown in the previous 
paragraph. With = 1.37 mm and Ve-n = 
1.92 * 10^ l/s we obtain De = 36.1 m^/s. 

Furthermore fi is related to by the Einstein- 
Smoluchowski equation (Eq. [t]) [1 . 



ePe 
K T 



(7) 
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K the Boltzman constant and T the hot electron 
temperature. With T = 25 keV we obtain /i = 
1.44 * 10~^ m? /sV. We can now compare the two 
terms on the right side of Eq. |5j 



lineE = 5.28* 10^ 



De Vne = 4.39 * 10^ 



1 



(8) 



(9) 



E is the electric field assumed to be the plasma 
potential of 3 F divided by the plasma chamber ra- 
dius R = 12 mm. Vng has been estimated by the 
electron density Ue divided by the plasma chamber 
radius. From the obtained values we can see the 
diffusion process is dominated by the term De Vng 
rather than the term jirie E. This allows simplify- 
ing the calculation of the number of collisions with 
neutrals an electron undergoes before being lost by 
diffusion (Eq. [lO]). 



^coll 



Td ye 



(10) 



Td is the mean time the electrons spend in the 
confinement before being lost by the discussed dif- 
fusion process. Ve-n is the electron collision fre- 



quency with neutrals. We can write Eq. [TO] as Eq. 

m 



ncoii 



rie Vpi 



asma ^e — n 



flDe Apia 



(11) 



Vpiasma IS the plasma volume, fiDe the electron 
loss rate equal to the right hand side of Eq. [9] and 
Apiasma the plasma surface. We approximate the 
geometry of the ECR plasma by a cylinder with 
length and diameter of the ECR plasmoid. Further 
introducing cva = Vpi 

asma 

/Api 

asma ^s the ratio of 
the plasma volume and the plasma surface leads to 



Eq. 12 



^coll 



Dp Up r? 



(12) 



With our values one obtains ricoii = 19.1. This 
means hot electrons with a kinetic energy of 25 keV 
undergo about 19 collisions with neutrals before be- 
ing expelled from the plasma by diffusion across the 



magnetic field lines. Relating to the CO2 gas oper- 
ation discussed further down a similar calculation 
for ions yields an average of 22 collisions with neu- 
trals before a singly charged carbon ion is lost by 
diffusion. 

We conclude from this section that the calcu- 
lated mean free paths of the hot electrons are long 
enough to fit many times into the plasma container 
and that diffusion across the magnetic field lines is 
not impeding the trajectories of the hot electrons 
because many collisions are needed to complete a 
successful diffusion loss mechanism. Hence the tra- 
jectory integration based on the simulated mag- 
netic field presented in the next section is valid. 
Furthermore, considering diffusion across magnetic 
field lines, ions are approximately as well confined 
as electrons in terms of number of collisions before 
loss. 

3 The magnetic field 

For ah magnetic simulations of SWISSCASE, the 
software TOSCA by vectorfields [4], a finite element 
solver, is used. A summary of the essential param- 
eters of the finite element model are presented in 
Table 



number of linear elements 


4.7* 10^ 


number of quadratic elements 


4.8* 10^ 


number of nodes 


1.6* 10^ 


number of edges 


6.1*10^ 


RMS 


1.033 


distortion of worst element 


1.40 * 10-^ 



Table 3: Some characteristics of the finite ele- 
ment model used for the numerical calculation of 
the magnetic field. 

Figure [3] shows the finite element model used for 
the magnetic field simulation. Cut planes have 
been introduced for the visualization of section 
views. 

The axial magnetic field has been measured using 
a Hall probe. The simulation result did not devi- 
ate from the measurement more than 2.8 % at any 
point along the z — axis (see Figure |3|. The maxi- 
mal error near the center of the electron cyclotron 
resonance zone was 0.97 %. Figures |4j Figures [s] 
and [6] present section views of the magnetic model 
defined by the cut planes introduced in Figure [3j 
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Figure 3: All magnetic parts of the finite element 
model have been meshed carefully according to size 
and relevance thereby optimizing the simulation ef- 
ficiency. 

The ECR quahfied surface, cahed plasmoid, 
presents itself compact and point symmetrically de- 
formed as expected from the permanent magnet 
setup. Figure [4]3 shows, despite its larger size than 
the permanent magnetic ring to the right (positive 
z-coordinate), the permanent magnet ring to the 
left (negative z-coordinate) does not deform the iso 
contour line in any asymmetric way. However the 
difference in size of the larger permanent magnet 
ring shifts the center of the plasmoid along the z 
axis toward the larger permanent magnetic ring (in 
Figure [4]3 to the left) without distorting its overall 
shape on a larger scale. 

4 The electron distribution 

In addition to the solution of the magnetic model, 
electron trajectories have been calculated using 
the finite element model trajectory integrator |4| 
SCALA by vectorfields. The initial velocities of the 
particle trajectories were chosen to be anisotropi- 
cally distributed. The electron velocity vector com- 
ponents perpendicular to the local magnetic field 
lines are distributed with a temperature of 25 keV. 
The velocity vector components parallel to the lo- 
cal magnetic field are distributed with a tempera- 
ture of 2 eV [8], the temperature of the cold elec- 
trons. This anisotropic velocity distribution com- 
plies with the ECR heating mechanism with the 




(b) 

Figure 4: Simulated magnetic field in plane A. For 
the purpose of orientation the outlines of the per- 
manent magnet structure and the Armco soft iron 
plate are represented by white lines in Figure a. 
Figure (b) features a detailed view of the center 
section. The ECR iso contour line is highlighted, 
revealing the topography of the ECR qualified sur- 
face. 

same anisotropic orientation which is perpendicu- 
lar to the local magnetic field lines ^ . In total 
36328 trajectories have been launched for a simu- 
lation run. The presence of space charge caused by 
the plasma electrons has been modeled by intro- 
ducing a charged sphere with the same diameter of 
the plasmoid and a potential of —3 V [3 [3 [9]. 

The magnetic field distribution of the MEFISTO 
ECR ion source located at the University of Bern 
has been simulated by Bodendorfer at al. [4 and 
matches the observed surface coating triangles in 
MEFISTO. The simulated magnetic field distribu- 
tion of SWISSCASE does not resemble the respec- 
tive surface coating triangles in SWISSCASE to the 
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Figure 5: Simulated magnetic field in plane B. 

same extent. However, a closer look at the electron 
distribution instead of the magnetic field distribu- 
tion, suggests a strong link between electron distri- 
bution and surface coatings in SWISSCASE. Fig- 
ures [7^-d illustrate how the shape of the simulated 
electron current density distribution changes from 
a hexagon at the plasmoid center into a triangular 
shape at the tip of the microwave antenna, where 
a triangular-shaped surface coating has been ob- 
served (Figure [8| after the operation of the SWISS- 
CASE ECR ion source [4 . This transformation of 
the electron density maps represents the changing 
cross section of the ECR zone while the cut plane 
is moved from the plasmoid center toward the an- 
tenna. Note due to the geometry of the magnetic 
field, the plasmoid and the electron density dis- 
tribution are point symmetric with respect to the 
plasmoid center. 

In the next section we present the observed coat- 
ing pattern and a spectral analysis of it which shows 
the coating consists indeed of particles found in the 
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Figure 6: Simulated magnetic field in plane C. Fig- 
ure (h) features a detailed view of the center section 
of Figure (a). The ECR iso contour line appears 
circular. 

plasma. 



5 The observed surface coat- 
ing pattern 

Figure [S] shows a typical surface coating triangle re- 
sulting from the exposure directly or indirectly to 
the ECR plasma. In this particular case the ob- 
served surface coating triangle was at the tip of the 
microwave antenna (see Figure [T]) in the SWISS- 
CASE facility after the operation of a CO2 plasma. 
In some cases the observed pattern is a sputtered 
area rather than a surface coating. This may origi- 
nate from the difference in particle energies the sur- 
face is exposed to. The observed surface coatings 
have been exposed directly to the plasma without 
any kind of extraction or post acceleration of the 
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Figure 7: Electron current density maps at different positions along the optical axis. Note the increasing 
similarity with surface coating triangles (see Figure^ toward the antenna tip at 15 mm offset from the 
plasmoid center. The white ring represents the plasma chamber wall (diameter 24 mm). Dimensions in 
mm. 



plasma particles. In this case the particle energies 
may have been low enough to enable a deposition 
process of the neutralized carbon atoms to form the 
observed surface coating. On the other hand sput- 
tering shapes are observed after the extraction and 
post acceleration where the involved particles have 
energies many orders of magnitude larger than they 
had inside the plasma. In this case the particle en- 
ergy might be too large for a continuous deposition 
process and hence sputtering occurs. 

The antenna tip was chemically analyzed to de- 
termine the composition of the surface coating. 
The result is summarized in Table [4j 



Element 


Weight [%] 


Abundance [%] 


Carbon 


90.61 


93.04 


Oxygen 


7.36 


5.67 


Fluorine 


1.41 


0.91 


Neon 


0.62 


0.38 


Total 


100 


100 



Table 4: Composition of the surface coating in mass 
and abundance fraction. 

The dominance of carbon in the coated surface 
suggests a deposition process originating from the 
plasma because the plasma mainly consists of car- 
bon and oxygen neutrals and ions. Oxygen is found 
as a minor fraction only because it mainly recom- 
bines with other oxygen and carbon radicals and 



molecules to form molecular oxygen, ozone, carbon 
monoxide or carbon dioxide. All of these recombi- 
nation products are gaseous and therefore not likely 
candidates for a deposition process. FC40, a liquid 
composed of Fluorine and Carbon, N-(C3F7)3, is 
used as cooling fluid enveloping the plasma cham- 
ber. Fluorine combines well with metals. There- 
fore, the Fluorine traces in the spectral analysis 
may originate from cooling fluid diffusing in small 
amounts through the copper plasma chamber wall. 
Neon found in the spectral analysis is a remainder 
of previous Neon plasma operations in the same 
plasma chamber. 

6 Relation between simula- 
tion and observation 

In this section we explain the relation between the 
observed surface coating pattern and the simulated 
electron distribution. 

As discussed in Section |2] the electron mean free 
path is limited to 4.9 m by collisions with neutrals. 
The mean free path is therefore far longer than the 
ECR plasmoid. This in turn potentially allows ion- 
ization processes anywhere inside the plasma con- 
tainer. However the magnetic simulation and the 
numerical trajectory integration show that the elec- 
trons are confined to a limited space inside the con- 
fining field (see Figures[7^-d for the electron density 
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Figure 8: Photograph of the microwave antenna of 
SWISS CASE after operation (CO2, 100 W input 
power, 10.48 GHz), published by Bodendorfer at al. 
in l4 ^ dark triangular patch consisting mostly of 
carbon is clearly visible, which was not present be- 
fore the operation. The shape of the deposited coat- 
ing matches the electron density map, as explained 
in the text and in Figure^ The pattern is displaced 
from the optical axis due to a misalignment of the 
microwave antenna caused by the antenna's length 
and weight. The misalignment did not cause any 
measurable decrease in the ion beam performance. 
The antenna has a diameter of 12 mm and the side 
length of the triangle measures 4-2 mm. 

profile). 

To investigate the spatial distribution of the ions 
we make the approximation that the ion-neutral 
cross section equals to twice the value of the hot 
electron-ion scattering cross section (8.5*10~^^ m^) 
and a neutral density equal to the value presented 
in Section 2 of this article (4.805 * 10^^ 1/m^). 
This leads to a mean free path length for ions of 
0.8 m which is 62 times the length of the ECR 
plasmoid. This suggests the ions would completely 
fill all of the vessel volume resulting in an uniform 
density distribution. However in a magnetic field 
of 388.6 mT singly charged carbon ions with a ki- 
netic energy of 2 eV gyrate with a Larmor radius of 
1.82 mm around the magnetic field lines on which 
they were created. The trajectories of the carbon 
ions are therefore restricted by the magnetic con- 
finement field rather than the limited geometry of 
the plasma vessel. Hence the ions do indeed travel 



for a long mean free path length relative to the ECR 
plasmoid dimensions. But the ions are prevented 
by the magnetic field from filling all the vessel vol- 
ume and are confined to a restricted space similar 
to the electron distribution presented in the last 
section. 

This explains the good agreement between the 
observed surface coating triangle and the simulated 
electron distribution because the carbon ions do not 
travel far from the region of ionization defined by 
the presence of high energy electrons. 

7 Conclusion 

We simulated in detail the 3D magnetic field dis- 
tribution of the SWISSCASE ECR ion source. 
The simulation of high energy electron trajecto- 
ries based on the same magnetic model leads to 
shapes nearly identical to the observed ion surface 
coating triangle. This implies the ion distribution 
is related to the high energy electron distribution. 
This is justified by the magnetic field parameters, 
the ion mean free path lengths and similar Larmor 
radii and diffusion parameters of hot electrons and 
singly charged carbon ions in SWISSCASE. 

Furthermore we generalize that lightweight ions 
with Larmor radii and mean free path lengths sim- 
ilar to those of a hot electron population can be 
considered as spatially bound to the latter. This 
can significantly simplify the calculation of the ion 
distribution in ECR related applications such as 
ECR ion engines and ECR ion implanters. 
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